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Abstract 



By combining Pauli algebra with distribution theory, we give a very 
compact and conceptuaUy simple formulation of Huygens' principle in 
classical electrodynamics. The Stratton-Chu and Kottler-Franz repre- 
sentations of an electromagnetic field as surface integrals are derived 
with minimal effort and maximal clarity. They are then generalized by 
allowing the integration surfaces to move freely, so that charge distri- 
butions in arbitrary motion are represented. 



1 Electrodynamics with the PauU algebra 

The Pauli algebra is a generalization of vector analysis in R^. One defines 
the products of two vectors A, A' by 

AA' = A-A' + iAxA'. (1) 

Thus A A' consists of a scalar part and a (pseudo-) vector part, the latter 
represented by an imaginary axial vector, which is actually the oriented area 
spanned by A and A': 

{AA')s = A- A', {AA')^ = iAx A'. 

The two common bilinear expressions are thus united into a single complex 
entity. A general element of the algebra, here called a Pauli number, is 
represented by a complex scalar plus a complex vector 

A = Ao + A, Ao G C, A€C^ 
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and the product of two such elements is 

AA' = AoA'o + A- A' + AoA' + AA'Q + iAx A' (2) 

A concrete representation of the algebra is given in terms of 2 x 2 matrices 
by the correspondence^ 



■^0 + As Ai + iA2 



^0,^1,^2,^3 e C, (3) 



where A A' is represented by the matrix product. Define the spacetime 
differential operators 

D = at + V and D = 5* - V, 

which act on a Pauli-valued field A(a;) = A{x, t) by 

BA={dt + V)(Ao + A) = (io + V • A) + (A + VAo + fV x A) (4) 
DA = {dt - V){Aq + A) = (io - V • A) + (A - VAq - iV x A), 

where Aq = dtAQ, A = dtA, and whose product is the scalar wave operator: 

DB = BP = - = □. 

Now consider the scalar wave equation 

Uf{x) = g{x) (5) 

where g{x) is a given source function which, for convenience, is assumed to 
be a distribution of compact support. The wave radiated by g is the unique 
causal ^ solution 

fix) = I d^x' P{x - x')g{x') = P * g{x), (6) 



^Along with complex numbers and quaternions, the Pauh algebra is one of the sim- 
plest examples of Clifford algebra. Although its first application to physics was in quantum 
mechanics, it has also turned out to bo useful in other fields, especially classical electrody- 
namics [B99]. The relation to Pauli matrices is found by choosing Ao = and A^: = Vxk 
to be the unit vector in the direction of Xk- The matrix representing Ak is then the Pauli 
matrix ak- 

^In this context, causality simply means that / is supported in the future region of 
g. If g vanishes aX t = — oo as assumed here, it sufiices to take the 'initial condition' 
f{x, -oo) = 0. 
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where * denotes spacetime convolution and P is the retarded propagator, 
which is the wave radiated by g{x) = S{x) = S{x)d{t): 

P{x) = P{x,t) = ^-^^^, aP{x) = 5{x). (7) 

Hence the wave operator is invertible on the space of such fields, with 

□"^ = P* . 

Since □ is a scalar operator, it operates on Pauli fields A(x), expressed in 
Cartesian coordinates, by 

DA{x) = DAo{x) + DA{x). 

Thus we may extend the wave equation (5) to Pauli fields as 

D¥{x)=G{x) (8) 

where G{x) is a Pauli- valued distribution with compact support. The unique 
causal solution is 



¥{x) = P* G{x) = [ d^x' P{x - x') G{x') 



(9) 



We now apply the Pauli algebra to classical electrodynamics, more or less 
following [B99]. To minimize the appearance of unnecessary parameters, we 
use natural Lorentz-Heaviside units, where Eq = fj,o = c = 1. An electromag- 
netic field in free space consists of two vector fields E{x),H{x) satisfying 
Maxwell's equations 

E -V X H = -J V ■E = p (10) 

H + VxE = V-H = (11) 

where {p, J) is a given charge-current density. The obvious symmetry of 
these equations suggest combining the two fields into a single complex field 

F{x) = E{x) +iH{x), (12) 

for which Maxwell's equations reduce to 

F + iVxF = -J V-F = p. (13) 
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Now interpret F{x) as a Pauli field with vanishing scalar component. Then 
(4) shows that (13) further reduces to the single equation 

lDF = p-J=S. (14) 

The homogeneous equations (11) state that the source J is real, but it will 
be useful to allow J to be complex: 

J = Jg + i ~ Pe ~ J e Jm ~ Pm ~ J m (l^) 

where Jg and represent electric and magnetic sources, respectively. Al- 
though Maxwell's equations require J„j = 0, virtual magnetic sources will 
be needed in the formulation of the general Huygens principle. 

To solve (14) for F, apply B: 

□F = BDF = BJ (16) 

and note that 

B J = (5t - V)(p - J) = (p + V • J) + (iV X J - Vp - j). (17) 

Since the left side of (16) is a pure vector field, the scalar component of the 
right side of (17) must vanish. This gives the continuity equation 

(DJ), = p + V- J = 0, (18) 

whose real and imaginary parts state that electric and magnetic charge are 
conserved. Assuming the initial condition F{x, — oo) = 0, we obtain the 
unique causal solution 

F = P*{BI)=B{P* I), (19) 
where the last equality follows because □ commutes with D and D. Thus 

F{x) = BA(x) (20) 

where the Pauli field 

A = P*J = $- A with $ = P * p and A = P*J, (21) 
representing the 4-potential, is the causal solution of the wave equation 

□ A= B]D)A= BF = J. (22) 
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In fact, 



F= BA = ($ + V-A)- V$-A + iVx A (23) 

shows that A satisfies the Lorenz gauge condition^ 

# + V • A = 0. (24) 

If J is complex as in (15), then so is A: 

A = Ae + z A^ Ae = $e - Ae A^ = - A„. (25) 

The free Maxwell field is then given by 

E= Re F = -V$e - Ae - V X Am (26) 
ii- = Im F = - - A^ + V X Ae. 

Of course, the homogeneous Maxwell equations require = A^ = 0. But 
the expressions (26) with a virtual magnetic 4-potential A^ will be used to 
formulate Huygens' principle. 

2 Electromagnetic Huygens principle 

The assumption that J(x) is compactly supported was made for convenience 
and can be relaxed. While it is reasonable to assume that the spatial support 
of J is bounded at any time, we want to allow sources persisting in time, for 
example a set of charged particles following world lines or extended charged 
systems evolving in time. This includes, among other things, time-harmonic 
systems. The above results remain valid provided the integrals converge. 

Let the sources be spatially bounded. To simplify the analysis, assume that 
the spatial support of J(£C, t) is contained in the interior of a closed surface 

5 C at all times t.^ We assume that S" is a smooth manifold, at least of 
class C^. Denote the exterior of 5" by F and its interior by E'. Let \{x) be 

^Evidently, the Lorenz gauge is selected by causality. 

■^This will be generalized to sources in arbitrary motion in Section 3 by allowing A to 
depend on time. Here we assume a fixed surface S, as is commony done in the derivation 
of Huygens' principle. 
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a function such that^ 

xeE ^ X{x) > (27) 
X e S =^ A = and |VA| = 1 
xeE' ^ X{x) < 0. 

The characteristic functions x and x' of ^ ^' may be written in terms 
of the Heaviside step function if as ^ 

X{x) = HiXix)) = I, X'ix) = H{-Xix)) - 

\o, X e E' 

Then 

N{x) = Vx{x) = -Vx'ix) = 6s{x)n{x) (28) 

where 

n{x) = VX{x) and Ss{x) = H'iXix)) = S{X{x)). 

Thus n is the outward unit normal on S and d^xds{x) is the 2D area 
measure on S: 

d^xdsix) = dS{x). (29) 

Remark: Since the characteristic function x{x) does not depend on the 

choice of A, neither does the distributional field N = Vx- The introduction 
of A is thus seen to be merely a convenience which makes the concepts easier 
to understand by using the relation H' = S. Similar remarks apply when 
X{x,t) is time-dependent, allowing for moving boundaries. 

Let F' be an interior field whose source 

/ = BF' (30) 



®An example of a function with these properties is 

{d{x), X £ E 
0, a; e S 
-d{x), X € E'. 

where d(x) is the shortest distance from x to S. 

®For a; € iS, we define x{x) = 'x!{x) = 1/2; but this singular case will not be needed in 
the sequel. 
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is spatially supported in the exterior region E. (This includes the case 
3' = 0, where F' is globally sourccless.) Since the spatial support of J' is by 
definition closed, it must actually be contained in some closed set V C E. 
Hence F' is defined and sourceless in an open neighborhood of S as well 
as in its interior E'. Thus both F and F' are defined and sourceless on a 
neighborhood of S. 

We shall construct a field F'^ whose soTirccs arc concentrated on S at all 
times and which coincides with the given field F in E and with F' in E'. 
The two partial fields arc 'glued' into a single field defined by 

F%x) = xix)F{x) + x'ix)F'{x), X = (x, t) € M^, (31) 

and the source of F^ is defined by applying B in a distributional sense: 

= p'-J'= BF^. (32) 

Note that Bx = Vx = AT and 

B(xF) = ( Bx)F + X BF = NF + x OF. 

Similarly, since Bx' = Vx' = — Vx = —N, 

B(x'-F') = -NF + x'^F'. 

Therefore 

F = JVF^' + XJ + X'/ (33) 

where 

= F-F' = + iW = E- E', = H - H' (34) 

is the jump field across S. Since J is supported in E' and J' is supported 
in E, we have the global identities 

XJ = and x'J' = 0. 

Hence 

js ^ = N.F^ + iN X F^ = Ss{n ■ F^ + in x F^) (35) 

is a distributional (4D) charge-current density supported on S, with a surface 
charge-current density {a, K) given by 

= 5sa a = n-F^ (36) 

J^ = dsK K = -inxFK 
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Like J, J'^ satisfies the distributional continuity equation 

(BJ^), = /9^ + V- J^' = 0, (37) 

which states that charge, now restricted to flow on S, is conserved. 

Note that even though JI is real, is in general complex, consisting of 
electric and magnetic sources on S: 

J^=Jf+iJ^ Jf=pf-Jf rm = Pm-Jm (38) 

with 

pt = dsae, ae = n-E^ Ji = SsK,, K, = nxW (39) 

Pm = ^S(^m, (Tm = n- W = 5sKjn, Km = -U X (40) 

If we wish to construct a physically realizable surface source J"^, then the 
absence of magnetic monopoles requires it to be real: 

J;^ = O n-W=Q and n x = Q on S. (41) 

That is, the normal component of and tangential components of E^ 
must be continuous across S. It can be shown^ that the scalar condition 
follows from the vector condition and Maxwell's homogeneous vector equa- 
tion. Since we are free to choose any sourceless interior field F' , (39) and 
(41) can be viewed as a set of boundary conditions for {E' , H') with {E, H) 
given. Thus we look for an interior field F' = E' + iH' such that 

f' + iV X F' = in E' 
nx E' = nxE and n H' = n H on S. (42) 

(Recall that F' actually extends as a sourceless field to a neighborhood of 
S.) This boundary- value problem has a unique solution if F is continuous 
in an open neighborhood of S, which will be the case if J is continuous in 
time.^ (Recall that we have also assumed S to be of class C^.) For this 
unique interior field, (39) and (41) are the jump conditions on the interface 
between the interior and exterior regions [J99, pp 16-18]. 

If the interior field does not satisfy (42), then the Huygens representations 
we are developing, while useful mathematically for expressing the given 'real' 

'^This follows in the frequency domain from Equation (6.38) in [CK92]. 

*This is a sufficient but not necessary condition, as follows from the properties of the 
propagator (7). Due to the factor S{t — r), the spread of J in both time and space tends 
to smooth F. 



8 



field F {i.e., with = 0) in terms of surface integrals, cannot be realized 
physically by actual surface sources. This is what was meant by saying that 
the magnetic sources are 'virtual.' In either case, we now derive the 

Huygcns representations. 

The 4-vector potential for F^ is given by (21) and (29) as 

<t^{x) = P*p^{x) = [ d'^x' Ss{x')P{x - x')n{x') ■ F^ {x') (43) 

A^{x) = P*J^{x) = -if d^x' Ss{x')P{x - x')n{x') x F^{x') 

where r = \x — x'\ and 

[u] = u{x', t — r), [v] = v{x', t — r) 
denote retarded expressions. Hence the electric and magnetic 4-potentials 



are 



$f (x) = ^J dSix') r-^ [n ■ E^] (44) 
Al{x) = ^j dS{x')T-^ [nxW] 
^U^) = ^ J dSix')r-' [n-W] 
A'^{x) = -^ I dS{x')r-' [nxE^. 



Substituting these into (26) gives 



(45) 



If we choose F' = and assume x e E, then F^ = F and (45) reduce to the 
Stratton-Chu equations [HY99, page 32]. Since F' = does not satisfy (42) 
except in the trivial case F = 0, the Stratton-Chu formulation of Huygens' 
principle requires virtual magnetic sources. 

However, if F' is chosen to be the unique solution of (42), the Stratton-Chu 
equations reduce to the simpler expressions 

E^ = -V$f -Al = V X Af . (46) 
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Returning to the general case (44) and (45), note that involves only 
the tangential components of on S while involves only the normal 
components. The latter can be eliminated as follows. Begin with 

= -iV xF^ -J^ = -iV X {iV X - - A^) - 
= V xV X A^ + iV X A^ + iN X F^. 

This involves only the tangential component n x F^ of F^ on S, and it can 
be integrated using the initial condition A^{x, — oo) = to obtain 

= V X V X 5^-^ + iV X A^ + iN X d^^F^ (47) 
where di^A%x,t)= f dt' A%x,t'). 

J —oo 

If we choose F' = and assume x & E, then (47) reduce to the Kottler- 
Franz equations [HY99, page 34]: 

= V X V X df^A^ - VxA^ l/^ = VxVx d^^A^ + V x . 

(48) 

However, note that (47) arc global, remaining valid when x ^ E' (where 
ps — pi"^ Sijid, in a distributional sense, even when a; G 6" as indicated by 
the last term. 

Like the Stratton-Chu equations, (47) involve virtual magnetic sources on S. 
If we assume that the interior field satisfies the physical boundary conditions 
(42), then (48) simplify to 

= V X V X dt^Al = V X Af . (49) 

3 Moving sources 

The above can be generalized to sources in arbitrary motion simply by letting 
A depend on time, so that the 2D surface 

St = {x: X{x,t) = 0} C 

enclosing the source S{x, t) at time t is time-dependent. The 3D hypersur- 
face 

S = {{x,t) : X{x,t) = 0} CM."^ 
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now represents the history of Sf. It is the oriented boundary separating the 
exterior and interior spacetime regions: 

E = {x: X{x) > 0}, E' = {x: X{x) < 0}, S = dE' = -OE. 

Define the Huygens field 

F^{x) = x{x)F{x) + x'{x)F'{x) (50) 

where 

X{x) = H{\{x)) and x'{x) = H{-X{x)) = 1 - x{x) 

are the characteristic functions of E and E'. Then the same arguments as 
above give 

js = ^ps ^ ( Bx)F^ = Ss{XF^ +n-F^ + inx F^) (51) 
with charge- and current distributions 

= 5sn- F^ , = -Ss {XF^ + inx F^) (52) 

where 

5s{x) = 6{X{x)) (53) 

is a distribution supported on S, whose interpretation is given by the mea- 
sure 

d^x5-s{x) = dtd?x5{X{x,t)) = dtdSt{x). 
The term —XF^ is a drag current on generated by its motion. Since 

-XF^ = An X (n X F^) - Xn{n ■ F^), 
it has both normal and tangential component. Equations (43) generalize to 

^%x)= f d'^x'ds{x')P{x-x')n{x')-F^{x') (54) 

A^{x) = - [ dV(5s(x')P(x - x') \x{x')F^{x') + in{x') x F^(x')| 

The electric and magnetic 4-potentials are the real and imaginary parts, 
and the Stratton-Chu equations for a moving surface are obtained exactly 
as in (45). Enforcing the homogeneous Maxwell equations on S gives the 
boundary conditions 

n-W=0 and XH^ +nxE^ = 0. (55) 

Note that the vector condition implies the scalar condition if A ^ 0. (As 
stated earlier, this is also true when A = 0, though not as obviously.) It can 
then be used to determine the interior field. 
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4 Local partitions and nonlinearity 



The expression (50) defines a global field using a partition of spacetime 
into the exterior and interior regions separated by the interface S: 

R'^ = EUSUE'. 

This can be generalized to a finite (or even infinite) sum of partial fields 
defined on regions with characteristic functions Xk- 

F{x) = ^Xk{x)Fk{x), (56) 



where the superscript S has been dropped. We assume that F^ is a source- 
less field in an open spacetime region Ok containing the closure of Ek, thus 
extending beyond its boundary. Since 

Xk^Fk = yk, 

the source of F is the distribution 

J = BF = ^( Dxfe)n = ^ {xkFk + Vxk ■ Fk + iVxk x Fk} (57) 
k k 

^ P = Y1 Vxit -Fk and J = {XkFk + i^Xk x Fk] . 

k k 

Since ^Xk is supported on dE^^ J is supported on the 'cellular' boundary^ 

S = [jdEk. (58) 
k 

Furthermore, since^'' 

xedEkndEi B{xk + Xl) = 0, (59) 
J depends only on the jump fields 

Fii = Fk-Fi 

® Since dEk is oriented by the unit normal pointing into its interior Ek, each interface 
dEh n dEi between adjoining regions occurs twice in (58), with opposite orientations. 
Hence the oriented sum (chain) dEk vanishes but the set-theoretic union S does not. 

'^"it makes no sense to say that Dxfc = — Dxi since Dxfc and Dx( are both infinite on 
dEkC\dEi. On the other hand, Xfc+Xi is the characteristic function of Ek^EiU{dEkr)dEi), 
hence (59) makes sense. 
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across the interfaces between adjoining regions Ek,Ei. Thus 

F=nA where A = P*J (60) 

gives a generahzed Huygens representation of F in terms of sources sup- 
ported on S. Furthermore, the projection property 

Xl{x)Xm{x) = SimXlix) (61) 

imphes that nonhnear expressions in F have similar partitions. For example, 
the scalar Lorentz invariant 

F'^ = F ■ F = E"^ - H'^ + 2iE ■ H 

has the local partition 

k 

and the electromagnetic energy-momentum density 

S = ^FF* = ^{F ■ F* + iF X F*} = U + S (62) 
U = ^{E^ + H^), S = ExH 

has the local partition 

k k k 

Hence the local power density (rate of increase of energy density) is 

(B§), = f/ + V- S = ^ iXkUk + Vxfc • Sk} + ^Xfc {f/fc + V • Sk] . 

; k 

Since F^ is sourceless in Ok, it follows from Poynting's theorem that 

Uk + V-Sk = in Ok 



Equation (61) fails numerically on dEkCidEi where Xkix) = Xi{x) = 1/2, but it holds 
weakly, in the sense of distributions, i.e., 

j A'^xxi{x)xm{,x)f{x) = 6im j d'^x Xi{x)f{x) 

for any continuous function / with compact support or rapid decay (needed when Ek or 
El are unbounded). 
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and thus 



if + v ■ s = Y,{xkUk + Vxk - Sk} . 



(63) 



Here XkUk is the rate of increase in the energy density coming into due 
the motion of the boundary dE^, and Vx^ • Sk is that due to the incoming 
momentum flowing through dEk. Due to (59), the right side of (63) involves 
only the differences 

Ui, = Uk-Ui and Si, = Sk-Si on dEkHdEi. 



Since the general partition (56) allows arbitrary choices of sourceless fields 
Fk in domains Ok containing the closure of Ek, these differences need not 
vanish. By enforcing boundary conditions on any interface dE^ fi dEi, the 
corresponding terms can be made to vanish. But then that interface can be 
removed, thus merging the two cells into one. 

It is instructive to confirm (63) using the expression (57) for the surface 
current J. Recall that J is generally complex, including magnetic as well as 
electric sources: 

The generalized Poynting theorem for a complex surface current density is 
derived by applying the distributional Maxwell equations 



F + iV X F = - J 



F -iV xF* = -J* 



to 



U + V . S =^^F ■ F* +F ■ F* +iV X F ■ F* -iF -V X F*^ , 

which gives 

U + V-S = ~{J-F* + F-J*) = -Je-E-Jm-H. (64) 

The right side is, like J, a distribution supported on S. The partitions 
- J = ^ {XfeFfc + iVxfc X Ffc} F = ^xiFi 
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give 

-J F* = Y, {XlXkFk ■ F*i + ixiVxk X Fk ■ F^} . 

kl 

Using Vxk X Fk ■ F^ = Vxk ■ Fu x Ff, (64) gives 

U + V-S=]^Y. i(x«Xfe + XkXl)Fk ■ F*i + iixNXk + Xk^Xl) ■ Fk x F\} . 
kl 

But the projection property (61) implies the distributional identities 

XlXk + XkXl = hlXk and Xl^Xk + Xk'^Xl = hNXk, 

therefore 

U + V-S = Y, {XkUk + VXk ■ Sk} 

k 

in agreement with (63). This confirms the consistency of our computations 
involving bilinear distributional expressions.^^ It gives reason to hope that 
partition methods could be useful in the synthesis of global solutions to 
nonlinear equations from local solutions. 
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